Quantifying Nonequilibrium Behavior with Varying Cooling Rates 
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We investigate nonequilibrium behavior in (l+l)-dimensional stochastic field theories in the context 
of Ginzburg-Landau models at varying cooling rates. We argue that a reliable measure of the 
departure from thermal equilibrium can be obtained from the absolute value of the rate of change 
of the momentum-integrated structure function, AStot ■ We show that the peak of AStot scales with 
the cooling, or quench, time-scale, r q , in agreement with the prediction by Laguna and Zurek for 
the scaling of freeze-out time in both over and under-damped regimes. Furthermore, we show that 
the amplitude of the peak scales as r 9 ~ 6/ ' 5 independent of the viscosity. 
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Due to the widespread role of symmetry break- 
ing phase transitions in cosmology, condensed matter 
physics, and high energy physics, it is imperative that the 
nonequilibrium aspects of the dynamics of these phase 
transitions be better understood. In particular, the effect 
of the cooling rate is still a widely open topic of inves- 
tigation. Since all but the simplest problems are noto- 
riously difficult to handle analytically, most approaches 
assume one of two extreme limits: either one assumes 
a quasi-adiabatic (infinitely slow) cooling, or an instan- 
taneous (infinitely fast) quench, often even in numerical 
studies. However, in reality all cooling occurs with a fi- 
nite rate and, therefore, lies somewhere between these 
two extremes. Some notable exceptions are Wong and 
Knobler's 0] and Binder's || work on quenches within 
the one-phase region, Wong and Knobler's || and Ruiz's 
j| work on double quenches, and Onuki's work on 
periodic quenches. For excellent reviews, see ||,jj7J and 
More recently, Zurek, Laguna and Zurek (9jT and 
Bettencourt et al. jl(| have investigated the effect differ- 
ent cooling rates have on defect density in classical fields, 
while Bowick and Momcn have examined this issue 
for quantum fields. 

Although these works have done much to advance our 
understanding of this issue, the effect of different cool- 
ing rates on the dynamics of phase transitions, as well 
as other nonequilibrium situations that commonly arise 
in nonlinear field theories, requires much further study. 
Even if advances in computer technology in recent years 
have greatly facilitated numerical studies of phase tran- 
sitions, some aspects of nonequilibrium dynamics of non- 
linear field theories remain difficult to study numerically, 
as they necessarily require large lattices and/or long sim- 
ulations. In particular, dynamical studies near the criti- 
cal region of continuous phase transitions, or of nuclc- 
ation in discontinuous phase transitions, require both 
long runs and large lattices. 

Once we focus on dynamical issues, it is often desirable 
to observe how a system evolves towards its final equilib- 
rium state. It is thus important to develop tools designed 
to quantify the nonequilibrium behavior of a system be- 



ing cooled (or warmed up) in a way that is numerically 
efficient. With this goal in mind, in this letter we propose 
a possible measure to quantify the departure from equi- 
librium of a system coupled to a heat bath, which clearly 
correlates the approach to equilibrium with the relevant 
control parameters of the system, namely, the absolute 
value of the rate of change of the momentum-integrated 
structure function, A»S tot , which will be defined below. 

For simplicity of comparison, we chose to use Laguna 
and Zurek's model of a linear pressure quench for a single 
classical real scalar field <f>(x,t) in (l+l)-dimensions [p|. 
Although there is no phase transition in one dimension, 
we can still examine (local) symmetry breaking and phase 
ordering through the formation of kink-antikink pairs. 
The present work should be considered a further step in 
the quantification of nonequilibrium behavior, which can 
be extended to higher-dimensional systems. All quanti- 
ties have been scaled appropriately to be dimcnsionless. 

We implement linear cooling with the following poten- 
tial (or, equivalently, the homogeneous part of the free- 
energy density for the order parameter 6), 
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V{6) = - (1 - 2e(t)6 2 + 6 4 ) 



(1) 



where e(t) = min(i/r 9 , 1), and t is measured from the 
"phase transition point" (when V"(6 = 0) = 0), so that, 
when the potential is a single well, both t and e(t) are neg- 
ative. r q is the cooling time-scale. Note that the poten- 
tial stops changing when e(t) = 1, which, in a Ginzburg- 
Landau system, is the zero-temperature limit. We simu- 
late the coupling of the scalar field 6 to the thermal bath 
using a generalized Langevin equation, 
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The terminology "pressure quench" is justified by the 
fact that the temperature of the bath T, which is related 
to the viscosity rj and the stochastic force of zero mean 
£,(x,t) by the fluctuation-dissipation relation, 



(£(x, t)£(x', t')) = 2i]TS{x - x')5(t - t') 



(3) 
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remains constant, while the quadratic coefficient of the 
potential changes linearly in time. This is equivalent to 
a pressure quench at constant temperature in the labora- 
tory. An obvious extension of this work is to implement 
a true linear cooling, where the temperature of the bath 
changes linearly, and e(t) = (T c — T)/T c . We will leave 
this case for a future investigation. [See, however, Ref. 
JTo| .] Our main interest here is in proposing a measure for 
nonequilibration of classical fields, which can be adapted 
for several different situations, including those involving 
cooling through the bath temperature, as we show below. 

Zurek, and Laguna and Zurek have studied the effect 
of different cooling rates on the density of zero crossings 
of the field. [A counting of the number of times the field 
goes through zero for a given lattice length), which pro- 
vides the approximate kink density as a function of time 
H .] In general, the field will attempt to keep up with the 
changing potential as best it can, that is, its modes will 
try to keep thermalized as the cooling occurs. Clearly, as 
the cooling rate is increased, we can envisage a situation 
where this will not be possible any longer, and the field 
becomes "frozen" , unable to maintain thermal equilib- 
rium with the bath. According to Zurek's conjecture, this 
freeze-out occurs when the dynamical relaxation time - 
given by ~ \4>/4>\ for overdamped systems, and by 

~ I0/0I 1 / 2 for underdamped systems - is comparable 
to the time to (from) the phase transition. From this 
freeze-out condition, one can derive the scaling relation- 
ships for the freeze-out time, oc Tq 1 / 3 and oc T q x / 2 . 
Using these results, Zurek and Laguna find scaling laws 
for kink density, which they confirm with simulations and 
contrast with experimental results for pressure quenches. 
Here, we are mostly concerned with how to extend our 
knowledge of nonequilibrium properties of field theories. 
We will thus be using Zurek and Laguna's model as a 
testing ground for our methods, comparing some of our 
results to theirs. 

For any given moment during the evolution of the sys- 
tem, there will be local fluctuations around the space- 
averaged order parameter. These can be studied with 
the structure function Sk(t), which tells us how different 
Fourier modes evolve in time. Its time derivative will thus 
give us information on the rate of change of the individ- 
ual modes. (Notice that this is not the dynamic structure 
function, which is a Fourier transform in both space and 
time.) Integrating the time derivative of the structure 
function over wave number gives the net change of the 
fluctuations in the order parameter. For simplicity, we 
use the absolute value of the momentum-integrated time 
derivative of S&(i); as we will soon see, it will display 
a peak which will give us valuable information on the 
nonequilibrium dynamics of the system. We thus define 
the quantity AStot as , 
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Sk (t) — \uk (t)\ 2 , where u k (t) is the Fourier transform of 
the field fluctuations, u (x, t) = <f) (x, t) — (j) (t) and 4> (t) is 
the spatial average of the field. Sx is the lattice spacing 
and L is the lattice size. The integral is over the first 
Brillouin zone. Here is our main point: In equilibrium the 
field is in a steady state, and so AS to t ~ 0. Thus, values 
of AS'tot greater than zero can be used as a measure of 
nonequilibration of the system. AS'tot may be smoothed, 
if necessary, by sampling every few time steps, which is 
equivalent to averaging over the same number of time 
steps. 




V(0)=l/2 s 

AS lol (sim.,T„=.01) 
_ 4S M (theo.,T =.01) 
. _ KED (sim.,T =.01) 
AS lol (sim.,T =. 005) 
iS^ (theo.,T =. 005) 



FIG. 1. AStot and the average kinetic energy for an instan- 
taneous temperature quench in a time-independent potential. 

As a first illustration, Figure [I] shows AS to t and the 
kinetic energy density, KED, for a field in a time- 
independent potential, V(cf>) — \4> 2 ■ 4> is initially ther- 
malized to a temperature To with viscosity 77, which is 
equal to unity for both cases shown. At t = 0, the bath 
is instantaneously quenched to zero temperature. No- 
tice how initially AStot is approximately zero because 
the field is thermalized (KED = To/2), but begins im- 
mediately to rise when the quench begins, reaches a 
maximum and returns to zero after the field thermal- 
izes (KED — * 0). This case can be solved analytically 
for 77 < 2, if we take the ideal limits Sx — > and L — > 00, 
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where uj = y 1 — Jo,Ji and J 2 are Bessel functions 
of the first kind and iF% is a hypergeometric function. N 
and Z? = -^r, where C is the value of dS g^ at t = 0, are 
parameters that were varied to fit the data. This is shown 
for both initial temperatures, To = 0.01 and To = 0.005, 
as a line in Figure 0. The simulation data were fit using 
N = 0.82 and D = 0.65 for both T = 0.01 and T = 
0.005. We have verified (not shown) that both the time, 
^pcak, and amplitude of the peak, A poa ^ depend on the 
viscosity. Also, and this is very important, the amplitude 
increases linearly with temperature change, as can be 
seen in Figure [l]; thus, the location of the peak gives 
a measure of the equilibration time-scale of the system, 
while its amplitude provides a measure of the departure 
from equilibrium. 
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FIG. 2. AStot f° r different viscosities and cooling 
time-scales with the potential in Eq. 1. j4 pe ak increases with 
increasing cooling time-scales and decreases with increasing 
viscosity. £ pC ak/T"g increases with decreasing cooling time-scale 
and viscosity. 

In addition to testing our second order staggered 
leapfrog Langevin code for stability with the same pa- 
rameters as j| (L = 2048, Sx = .125, 5t = .025 and T = 
0.01), we reproduced their defect density results to test 
our code. For a more thorough description of the al- 
gorithm used see We added the calculation of the 
structure function to our code using a FFT routine from 
Numerical Recipes and the change in the structure 
function using a simple finite difference method, with the 
ability to average over a few time steps to smooth the 
data. We can then extract the amplitude and time (after 



the phase transition) of the peak. Figure ^ shows the 
variation in AStot(V r g) with viscosity and cooling time- 
scale. t/r q is used on the time axis for easier comparison 
between different cooling time-scales. The amplitude of 
the peak is higher the lower the viscosity and the shorter 
the cooling time-scale (faster cooling). Also, the peak 
occurs earlier for longer cooling time-scales (slower cool- 
ing). Those peaks that occur during the cooling show 
a plateau beginning after the peak and lasting until the 
cooling ends at e = 1. We were able to identify a peak up 
to T q ~ 128 at which point the data becomes too noisy 
to distinguish the peak from the plateau. The fastest 
coolings (with rj = 1 it happens for r q < 16 ) actually 
peak after the potential stops changing (e = t/r q = 1): 
the system remains out of equilibrium during the whole 
cooling process, which is thus equivalent to an instanta- 
neous quench. For r q > 256, AS to t — 0, and thus the 
system remains thermalized during the cooling; this is 
equivalent to an adiabatic cooling. It is between these 
two regimes that the cooling most dramatically affects 
the dynamics of the system. 

The time and amplitude of the maximum in AStot are 
plotted in Figure^ and Figure^ respectively as a function 
of r q for several viscosities. From fig. 3 we note that the 
location of the peak scales with r q in agreeement to the 
scaling obtained in Ref. Q for the freeze-out time. Since 
this is the time after which the field is able to begin 
relaxing to its equilibrium state, this seems reasonable. 
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FIG. 3. Peak time vs. cooling time-scale for 
time-dependent potential with varying viscosity. The straight 
lines above and below the data show the scaling predicted in 
Ref. [9] for over- and under-damped freeze-out times as a 
function of cooling time-scale. The simulation data interpo- 
lates between these two extremes appropriately. 
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FIG. 4. Amplitude of peak vs. cooling time-scale for 
time-dependent potential with varying viscosity. For all but 
the fastest quench times, all viscosities studied satisfy the 
scaling yl pca k oc t 9 ~ 6/ ' 5 . All three fit lines have a slope of -6/5; 
the fit lines are labeled by their y-intercepts. 

We found that the amplitude of the maximum univer- 
sally scales as t~ p , where p = 1.2 ± .05 for all but the 
fastest coolings. There are thus essentially three regimes: 
i) the field remains thermalized all through the cooling 
process, thus approximating an "adiabatic" cooling (for 
T q > 256 with r\ = 1); ii) the field remains out of equilib- 
rium through the complete cooling process, thus approxi- 
mating an instantaneous quench (for r q < 16 and r\ = 1); 
iii) the intermediate regime, where the field reaches ther- 
malization during the cooling process, which is signaled 
by the appearance of a plateau in AStot 7 which decreases 
in amplitude as r q is increased, until it cannot be dif- 
ferentiated from noise. The general picture can be de- 
scribed qualitatively as follows: the modes are initially 
thermalized, and their rate of change increases as they 
try to "catch up" with the changing potential (or, when 
appropriate, the changing environment). This may or 
may not happen during the cooling, determined by the 
value of T q . As the quench proceeds at the steady rate of 
r" 1 , AStot increases until the field begins to probe the 
bottom of the potential (the free-energy minima), when 
it has essentially "caught up" with the changing poten- 
tial. If this happens while the potential is still changing, 
the field then changes at the same rate as the poten- 
tial (AStot begins to plateau), until the potential stops 
changing and the field can fully thermalize (AStot ^ 0). 
In the fastest coolings (r q < 16 for rj = 1, the field is not 
able to "catch up" to the changing potential before the 
end of the cooling, and so AStot never plateaus; the am- 
plitude of the peak does not obey the scaling in Figure 



^ because the field is still changing after the cooling has 
ended, which, as we remarked above, approximates an 
instantaneous quench. We hope to report on the exten- 
sion of this method to higher dimensional systems and 
on a in-depth analysis of the relevant paramenter space 
in the near future. 
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